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Abstract 

(N " 

\ We study QCD at nonzero temperature and baryon density in the framework 

I of the analytic continuation from imaginary chemical potential. We carry out 

simulations of QCD with four flavor of staggered fermions, and reconstruct the 
phase diagram in the temperature-imaginary ^ plane. We consider ansatze 
for the analytic continuation of the critical line and other observables moti- 
vated both by theoretical considerations and mean field calculations in four 
CN \ fermion models and random matrix theory. We determine the critical line, 

' and the analytic continuation of the chiral condensate, up to /i^ — 500 MeV. 

■ The results are in qualitative agreement with the predictions of model field 

, theories, and consistent with a first order chiral transition. The correlation 

\ between the chiral transition and the deconfinement transition observed at 

\ fi = persists at nonzero density. 

I. INTRODUCTION 

r-| ! QCD at finite temperature and density is of fundamental importance, both on purely 

theoretical and phenomenological grounds. At high temperature asymptotic freedom will 
produce deconfinement and chiral symmetry restoration, at high density a richer phase 
structure and new phenomena have been predicted [1]. In principle, the lattice formulation 
provides a rigorous framework for the study of such phenomena. In practice, however, the 
lattice regularization is usually combined with importance sampling, which cannot be naively 
applied at nonzero baryon density, where the quark determinant becomes complex [2]. 

It has been recently realized [3] that this problem can be circumvented in the high T, low 
fi part of the QCD phase diagram where one can take advantage of physical fluctuations [4] 
Interesting physical information can be obtained by computing the derivatives with respect 
to fi at zero chemical potential and high temperature [5-7]. Fodor and Katz proposed an 
improved reweighting and applied it to the study of the four [8] and two plus one flavor 
model [9] . In refs. [10-12] the imaginary chemical potential approach was advocated and 
exploited in connection with the canonical formalism. In ref. [13] it was proposed that the 
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analytic continuation from imaginary chemical potential could be practical at high tempera- 
ture, and the idea was tested in the infinite coupling limit. In refs. [14] [15] the method was 
applied successfully to the dimensionally reduced model. In ref. [16] it was proposed that 
the critical fine itself can be analytically continued and results for two flavor of staggered 
fermions were presented. 

In this work we study QCD with four flavors of staggered fermions within the imaginary 
chemical potential approach. Some of the results presented here have been preliminarly 
reported in [17] 

In the next Section we review the formalism and the method. In Section III we recon- 
struct the phase diagram in the temperature-imaginary chemical potential plane. This is an 
interesting physical question by itself, it is a mandatory step toward the reconstruction of 
the phase diagram for real chemical potential, and provides some guidance for the analytic 
continuation. In Section IV we discuss a few aspects of the analytic continuation and offer 
two examples from model fleld theories: we will note there that by considering //^ rather 
than an analogy can be made between QCD at flnite baryon density and ordinary sta- 
tistical systems in external fields. The remaining part of the paper is devoted to numerical 
results at nonzero baryon density. The critical line is presented Section V, including a first 
assessment of the dependence on the number of flavors obtained combining the results by de 
Forcrand and Philipsen [16] with ours, and a cross check with the four flavor results by Fodor 
and Katz [9]. The results for the chiral condensate are presented in Section VI. In Section 
VII we discuss the nature of the chiral transition. Finally in Section VIII we summarize our 
results and give our conclusions. 



II. FORMALISM AND METHOD 



In the following we will briefly review the formulation of lattice QCD with a non-zero 
chemical potential /i and the possible uses of working with a purely imaginary /i. The zero 

e ~ j , with i^QCD the QCD Hamiltonian, 

can be discretized on an euclidean lattice with flnite temporal extent t — 1/T 

j {vUVi^Vi^) e-PSo[u\^-SAu,m ^ j (pu) e-f^^olu] ^gt M[U] (2.1) 

where U are the gauge link variables, ip and ip are the fermionic variables, Sq is the pure 
gauge action and Sp is the fermionic action which can be expressed in terms of the fermionic 
matrix M[U], Sp = ^M[[/]^. 

To describe QCD at flnite density the grand canonical partition function, Z{V, T, /i) — 

e T j , where N — J d^xip^ip is the quark number operator, can be used. The 

correct way to introduce a flnite chemical potential // on the lattice [2] is to modify the 
temporal links appearing in the integrand in Eq. (2.1) as follows: 

Ut — * e"^[/t (forward temporal link) 

U} e-'^'^Ul (backward temporal hnk) , (2.2) 
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where a is the lattice spacing. Sq is left invariant by this transformation but dctAf[?7] 
gets a complex phase which makes importance sampling, and therefore standard lattice 
MonteCarlo simulations, unfeasible. 

The situation is different when the chemical potential is purely imaginary: Ut — > e^°'^'Ut, 
Ut — > e~^"''^'Ut. This is like adding a constant U{1) background field to the original theory; 
det M[U] is again real and positive and simulations arc as easy as at = 0. The question 
then arises how simulations at imaginary chemical potential may be of any help to get 
physical insight in finite density QCD. 

One possibility is analytic continuation which should be practical at relatively high tem- 
perature [13]. Z{V, T, /i) is expected to be an analytical even function of /i away from phase 
transitions. For small enough /i one can write: 

log Z(fj,) = ao + 0-2/"^ + a^fj,^ + 0(^1") (2.3) 
log Z{fij) = ao - a2nj + a4nj + 0(/i^) . (2.4) 

Simulations at small hj will thus allow a determination of the expansion coefficients for the 
free energy and, analogously, for other physical quantities, which can be cross-checked with 
those obtained by reweighting techniques (see [18-20] for further material on the reweighting 
approach). This method is expected to be useful in the high temperature regime, where 
the first coefficients should be sensibly different from zero; moreover the region of interest 
for present experiments (RHIC, LHC) is that of high temperatures and small chemical 
potential, with fi/T ~ 0.1. This method has been already investigated in the strong coupling 
regime [13], in the dimensionally reduced 3-d QCD theory [14], and in full QCD with two 
flavours [16]. The Taylor expansion coefficients can also be measured as derivatives with 
respect to // at // = [5-7] . 

Z{V,T,ifij) can also be used to reconstruct the canonical partition function Z(y,T,n) 
at fixed quark number n [21], i.e. at fixed density: 



Z{V, T, n) = Tr ({e-^S{N - n)J = ^Tr (e'^ j^^ d^e^^^^"") 

1 /"^TT 



/ d^e-^^"Z(y, T, i^T) . (2.5) 

27r JO 

As n grows, the factor e~*^" oscillates more and more rapidly and the error in the numerical 
integration grows exponentially with n\ this makes the application of the method difficult 
especially at low temperatures where Z{V,T,ini) depends very weakly on /x/. The method 
has been apphed in QCD [11] and in the 2-d Hubbard model [10] [12], where Z{V,T,n) 
has been reconstructed up to n = 6 [12]. 

The study of the phase structure of QCD in the T - ifij plane is also interesting by its 
own, as we will discuss in the next Section, and to understand the ranges of applicability of 
analytic continuation. 

Results reported in the present paper refer to QCD with four degenerate staggered 
flavours of bare mass a • m — 0.05 on a 16^ x 4 lattice, where the phase transition is 
expected at a critical couphng Pc — 5.04 [22]. The algorithm used is the standard HMC 
algorithm. 
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III. THE PHASE DIAGRAM IN THE IMAGINARY ^ - TEMPERATURE SPACE 



Let us write Z{e) = Z{V,T,ieT) = Tr e^'^^e"^^ . Since N is a number operator, 



Z{6) is clearly periodic in 6 with period 2tt; moreover a period 27r/3 is expected in the 
confined phase, where only physical states with N multiple of 3 are present. However it has 
been shown by Roberge and Weiss [21] that Z{9) is always periodic 27r/3, for any physical 
temperature, and that the only difference between the low T and the high T phase should 
be a smooth, analytic periodic behaviour at low T, as predicted from a strong coupling 
calculation, and a non-analytic periodic behaviour at high T with discontinuities in the first 
derivatives of the free energy at 6 = 27r/3(A; + 1/2), as predicted from a weak coupling 
calculation. This suggests a very interesting scenario for the phase diagram of QCD in the 
T - i/ii plane which needs confirmation by lattice calculations. 

In order to get more insight into the phase structure of the theory, it is very useful to con- 
sider the phase of the trace of the Polyakov loop, P{x). Let us parametrize P{x) = \P{x)\e'^'^, 
and let (0) be the average value of the phase. In the pure gauge theory the average Polyakov 
loop is non zero only in the deconfined phase, where the center symmetry is spontaneously 
broken and (0) = 2A;7r/3, k — —1,0,1, i.e. the Polyakov loop effective potential is flat in 
the confined phase and develops three degenerate minima above the critical temperature. 
In presence of dynamical fermions P{x) enters explicitely the fermionic determinant and 
Z3 is broken: the effect of the determinant is therefore like that of an external magnetic 
field which aligns the Polyakov loop along (0) = 0. In the high temperature phase the Z^ 
degeneracy is lifted and (0) = is the true vacuum. 

When fij 7^ 0, what enters the fermionic determinant is P(x)e*^, 9 = fJ^i/T, instead of 
P{x). Therefore the determinant now tends to align (0) -|- 9 along zero, like a magnetic field 
pointing in the —9 direction. Hence one expects (0) = —9 at low temperatures. At high 
temperatures the fermionic determinant still lifts the Z^ degeneracy, but it is 9 that fixes 
which is the vacuum. In particular for [k — 1/2) < ^9 < [k + 1/2) one expects (0) ~ 2kn/3 
and 9 — 2{k + l/2)7r/3 should correspond to phase transitions from one Z^ sector to the 
other. 

In Fig. 1 we report our results for (0) versus the imaginary chemical potential for different 
values of (5. Since T = l/(iVja) and iV^ = 4 in our case, we have 9 = Aafii. For {3 = 4.94 
and 5.01, which are below the critical P ai /ij — 0, Pdl^i — 0) = Pc — 5.04, one has 
(0) —9 — —Aa/ii, i.e. (0) is driven continously by the fermionic determinant. For (3 — 5.10, 
which is well above f3c, we see that (0) ^ 0, almost independently of fij, as long as ^ < 7r/3, 
while for 9 > 7r/3 there is a sudden change to (0) ~ — 7r/3: we are clearly crossing the 
Roberge- Weiss (RW) phase transition from one Z3 sector to the other. At intermediate 
values, P — 5.065 and 5.085, (0) ~ until a critical value of a/ii, where it starts moving 
almost linearly with /ij crossing continuously the Z^ boundary: in this case there is no 
RW phase transition, but there is anyway a critical value of ni after which (0) is no more 
constrained to be ^ and changes again lienarly with 9: as we will soon clarify, this critical 
value of yU./ corresponds to the crossing of the chiral critical line, i. e. the continuation in the 
T-fij plane of the chiral phase transition. 

We display our results for the chiral condensate in Fig. 2. We expect a periodicity with 
period 27r/3 in terms of 9. Moreover (ipip), like the partition function, is an even function 
of iij: this, combined with the periodicity, leads to symmetry around all points 9 — mr/S, 
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with n an integer number, for {'ipip) as well as for the partition function itself. For f3 < /3c, 
('?/''?/') has a continuous dependence on a/// with the expected periodicity and symmetries. 
For P > Pc the correct periodicity and symmetries are still observed but the dependence is 
less trivial. At /3 = 5.065 there is a critical value a^i ~ 0.17 for which the theory has a 
transition to a spontaneously broken chiral symmetry phase: we are clearly going through 
the chiral critical line. The same happens for (3 = 5.085 at a/i/ ~ 0.22: in this case we 
have proceeded further, observing also the transition back to a chirally restored phase at 
ani ~ 0.30, which is, correctly, the symmetric point with respect to 6' = tt/S. At /3 = 5.10 
we never cross, when moving in /x/, the chiral critical line, but only RW critical lines ^. 

Besides the sets of runs at fixed /3 and variable we have also performed runs at fixed 
yU/ and variable f3 to look for other locations of the chiral line in the T-/i/ plane. In every case 
the location of the phase transition has not been determined by looking at susceptibilities, 
since our statistics for each single run were rather poor to this aim (of the order of 1000 
molecular dynamics time units) , but rather by looking at sharp changes of various physical 
quantities, among which the chiral condensate or the Polyakov loop. In each case sharp 
drops or jumps have been observed, allowing quite precise determinations of the transition 
point and suggesting the first order nature of the phase transition also at fij ^ 0. The drop 
of the condensate is always coincident with the sharp jump of the Polyakov loop, also at 
/i/ 7^ 0, suggesting that the coincidence of chiral symmetry restoration and deconfinement 
holds true also at /ij ^ 0. A summary of all our determinations of the chiral critical line is 
reported in Table I. 

It is interesting to illustrate in more details the determination of the endpoint of the RW 
critical line, f3E = 5.097(2). We have performed a simulation at exactly 6 = n/3, starting 
thermalization from a zero field configuration: we thus drive the system to one side of the 
RW critical line (assuming that the hue is there), i.e. on the border of one sector, since 
on the 16^ x 4 lattice it is already practically unfeasible to fiip through the RW line in 
reasonable simulation time. Let us now consider the baryon density, {b) = y^lnZ: it is 
an odd function of fi, since Z is an even function. Therefore, for an imaginary chemical 
potential, (b) is also purely imaginary and an odd function of /x/. This, combined with 
the periodicity in leads to the expectation that {b){9 = |~) = —{b){9 = I"*"). The 
last relation clearly implies that (6) = at = vr/S, unless (6) is not continuous on that 
point. Thus a non-zero value of (b) at 9 = implies the presence of the Roberge- Weiss 
critical line. On the right hand side of Fig. 3 the imaginary part of (6) at = |" is plotted 
as a function of (3: one can clearly see a transition from a zero to a non-zero expectation 
value, which permits the determination of Pe- We have verified that ai Pe also the chiral 
condensate and the polyakov loop have a sharp change, as we show in Fig. 4, and this implies 
that at this point, for = 7r/3, we also meet the chiral critical line, so that the RW critical 
line ends on the chiral critical line. On the left hand side of Fig. 3 we present instead the 
imaginary part of (b) as a function of /x/ for different values ot P < Pe- in this case (6) is 
always zero and continuous at ^ = |, but it is interesting to note how it starts developing 
the discontinuity as P ^ Pe- 



^ Error bars for the determinations at P = 5.10 and on the critical lines {6 = ir/3 and 9 = ir) are 
probably underestimated. 
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We present, in Fig. 5, a sketch of the phase diagram in the /^-/X/ plane, as emerges 
from our data and by exploiting the above mentioned symmetries. We can distinguish a 
region where chiral symmetry is spontaneously broken (indicated as IV in the figure) and 
three regions (I, II and III), which correspond to different sectors and repeat periodically, 
where chiral symmetry is restored. The chiral critical line separates region IV from other 
regions, while the RW critical lines separate regions I, II and III among themselves. As we 
have noticed above, the sharp drops fo the chiral condensate at the transition points suggest 
that the chiral critical line is first order also at /x/ 7^ 0, so that we expect all regions to be 
separeted by first order critical lines. 

IV. FROM IMAGINARY TO REAL AND THE CRITICAL LINE IN THE T, //^ 

PLANE 

We will concern ourselves with the properties of the critical line, as well with the ji 
dependence of a number of physical observables. 

To this end one needs to analytically continue the results from purely imaginary to real 
chemical potential: generically speaking, one deals with a function of a purely imaginary 
variable, continues it to the entire complex plane, and finally takes the limit e ^ of f{x+ie), 
X and e being real variables. General arguments guarantees that the analytical continuation 
is unique in the analyticity domain of the function. In practice, f{x) is unknown and has 
to be approximated by some series expansion or suitable ansatz. In either case, one has to 
pay attention to the fact that, by modifying the original expression by a non-leading term, 
the difference in physical quantities is still non leading. This can be complicated, and often 
mathematical arguments need to be supplemented by some physical insight [23], even in 
relatively simple cases [24]. 

The critical line itself can be analytically continued from imaginary to real chemical 
potential: a nice argument together with an application to the two fiavor model has been 
given in [16]. The analytical continuation of the critical line is easily discussed by considering 
that Z{fi), an even function of fi, is real valued for either real and purely imaginary yU : for 
real fi, the imaginary part of the determinant cancels out in the statistical ensemble, and 
it is even possible to cancel it exactly on a finite number of configurations by considering 
the appropriate symmetry transformation [25,26]; for imaginary chemical potential, the 
determinant itself is real, yielding authomatically a real partition function. For a complex 
H, Z{n) is, in general, complex. 

We can map the complex /x plane onto the complex /i^ plane, and consider Z{ii^) (on 
a finite lattice this is an exact polynomial). Then, Z is real valued on the real iJ? axis, 
complex elsewhere : the situation is analogous to e.g. the partition function as a function of 
a magnetic field, which becomes complex as soon as the external field becomes complex, and 
the physical domain (real partition function) is associated with real values of the couplings. 
The critical behavior of the system is then dictated by the zeroes of the partition function 
(Lee- Yang zeros) in the complex /^^ plane. The locus of the Lee Yang zeros is thought to be 
associated to a general surface of phase separation [27], and phase transition points, for each 
value of the temperature, are associated with the Lee- Yang edge building up in the infinite 
volume limit, thus defining a curve in the T, //^ plane. 
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This simple reasoning shows that it is very natural to think of the critical line as a smooth 
function T(^^), making it natural the analytic continuation from positive to negative /i^ 
values. Indeed, experience with statistical models shows that not only the critical line, but 
also the critical exponents are smooth functions of the couplings [28] (aside of course from 
endpoints, bifurcation points, etc.). Hence, they can be safely expanded, either via Taylor 
expansion or a suitable ansatz. In particular, /x^ = has no special character: it is just the 
point where the Lee Yang edge hits the real axis where T = Tc. 

Clearly the analytic continuation from imaginary chemical potential is practical when the 
critical hne is smooth, and a few coefficients suffice to describe it. It is then of some interest 
to discuss examples whose critical lines are exactly computable: we will indeed find that a 
second order expression in T and // well approximates the critical line over a large // interval. 

A. The Gross-Neveu model 

The Gross-Neveu model in three dimensions is interacting, renormalizable and can be 
chosen with the same global symmetries as those of QCD which, when spontaneously broken 
at strong coupling, produce Goldstone particles and dynamical mass generation. As such, 
the Gross-Neveu model (as well as other four fermion models) can provide some guidance 
to the understanding of the QCD critical behavior (see e.g.refs. [29-31]). 

The critical line for the three dimensional GN model was calculated in ref. [32] and reads 

1 - Af/Eo = 2r/Eo ln(l + e"'^/^) (4.1) 

where Eq is the order parameter in the normal phase. Setting = in the above equation 
gives the critical temperature at zero chemical potential, TcCa* = 0) = Eo/21n2 ~ .72Eo 
Expanding now ln(l + e^"^)) ~ ln2 — l/2x + l/8x^, and eliminating Eq in favour of 
we get 

(T - 1/27;)^ + /^7(8 In 2) = T'^/A (4.2) 

It is easy to check that this expression approximates very well the exact result 4.1, and 
then that a second order expression in is a good approximation to the critical line in this 
model. 

B. Random Matrix Theories 

As it is well known (see e.g. [33]), there is a remarkable relation between the symmetry 
breaking classes of QCD and the classification of chiral Random matrix Ensembles. 

For QCD with fermions in the complex representation (i.e. Nc > 2, fundamental 
fermions) with pattern of SSB SU{Nf)R x SU{Nf)L SU{Nf), the corresponding RMT is 
chiral unitary, /? = 2 in the Dyson representation. On the lattice, staggered fermions have 
unusual patterns of x^B: all real and pseudoreal representations are swapped. However, 
for complex representations, the corresponding RMT ensemble remains chiral unitary [33]. 
The critical line in the T, jj, plane for this ensemble derived in ref. [34] : 

(//^ + TY + - = (4.3) 
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is thus valid both on the lattice and in the continuum. Expanding it to 0{ii^) we obtain 



and by comparison with the exact result 4.3, we note that this expression describes well the 
critical line basically till its endpoint. 

V. THE CRITICAL LINE FOR THE FOUR FLAVOR MODEL 

In Section III we have presented our measurements of the critical points in the tem- 
perature imaginary chemical potential plane, see again Table I. Here we use those data to 
reconstruct the critical line for real chemical potential. We used both least ^^^^ least 
squares fits to a second order polynomial, and we studied the effect of a fourth order term. 

A significant sample of the results of the least ^^'^ least squares fits to a second order 
polynomials 

/?e(//) ^a + b/i^ (5.1) 

are collected in Table II and Table III. 

The weight of each data point for the least square fits is zero or one, i.e. we either include 
or discard it, and the quality of our fits -for different fi ranges indicated in the Table- is 
measured by the squared sum of residuals, to be compared with J^i^PY- Some of the fits 
(whose results we omit) included a linear term, and we confirm that it is compatible with 
zero, as it should on symmetry grounds. The least fits have been performed by discarding 
one point at the time. 

Prom the least square polynomial fits we obtain 

Pc = 5.038(2) + 0.94(7)/i^ (5.2) 

and from the minimum polynomial fits: 

Pc = 5.036(2)(2) + 0.98(3)(6)/x^ (5.3) 

The central values are the average of the results, the first error is the mean (absolute) 
deviation, the second one is the maximum error measured in individual fits. We combine all 
of these estimates to quote as our final result for the critical line in lattice units: 

(3 = 5.037(3) + 0.96(10)^1^ + 0(/x^) (5.4) 

The continuation to real chemical potential is shown in Fig. 6. The dotted fines are 
drawn in correspondence of the central values of the fit parameters 5.4 plus and minus the 
quoted errors. The dispersion remains reasonable in a fair interval of chemical potentials 
( /i = .2 corresponds, roughly, to a baryochemical potential of about 370 Mev and /x = .3 
to about 516 Mev.) Moreover, in the same interval, we note a nice agreement between our 
results and those obtained by Fodor and Katz via their improved reweighting [8] in the same 
four fiavor model and with the same quark mass = 0.05. Obviously our results, being 
uncorrelated, look less regular. 
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The line obtained via improved reweighting changes concavity around fi ~ .4. This is not 
easy to understand from a physical point of view, but, from a purely numerical perspective, 
this behavior could be reproduced by a negative fourth order coefficient in the expansion 
of the critical line. It would then be desirable to place at least bounds on the coefficient of 
the fourth order term, but, on the other hand, as the quadratic fits are satisfactory, this is 
not an easy task. To get a feeling of the effect of a fourth order term we constrained the 
constant and the second order term to leave the fourth order coefficient as a variable. The 
results are given in Fig. 7. The solid line corresponds to the result of a second order fit up to 
jjL — .2. The dotted fine is the correction to that result induced by a fourth order correction, 
which is very small. More sizable corrections are induced by constraining the second order 
coefficient to its extreme value (as inferred from the errorbars) and then fitting the fourth 
order term free parameter. 

From Fig. 6 and Fig. 7 we conclude that the imaginary chemical potential approach 
provides a safe estimate for the critical line up to ~ .3 (corresponding to a baryochemical 
potential of ~ 500 MeV) before loosing accuracy, and that within this interval imaginary 
chemical potential and reweighting give consistent results for the critical line. The loss 
of accuracy for the imaginary chemical potential results is mostly due to the influence of 
the fourth order term: we emphasize anyway that the results obtained with a fourth order 
term are much less accurate, however consistent with the ones coming from a second order 
polynomial. To reduce the error band we would need better data for the critical line at 
imaginary chemical potential. 

To convert to physical units we need the lattice spacing as a function of the coupling. 
We used as an input the lattice spacing measured dX (3 = 5.10 [35] to fix the scale in the 
two loop 13 function. We have verified that the ratio Vmeas of the scales taken at the two 
extrema of the interval of interest [5.04:5.10] measured in [35] agrees within a few percent 
with the ratio r^j coming from the two loop (3 function. Clearly the uncertainty induced by 
the interpolation and/or by the choice of the /? used as an input is less than the numerical 
errors on the lattice results. 

In Fig. 8 the external band show the results of the conversion to physical units by use 
of the two loop (3 function, the errors being those induced by the analytic continuation. In 
the same plot we draw the ellipse arc 



which turns out to be a nice approximation to the data in physical units, with k = 0.021 
(this just the result of the fit in the interval [0:400 MeV] to the central values of our critical 
line). Other functional forms would work as well, and, in particular, the parabola T = 
Tc{l — 0.5/c/i^/T^^) obviously provides a good approximation to the data, given the smallness 



The results by De Forcrand and Philipsen [16] combined with ours allow an assessment 
of the fiavor dependence of the critical line in QCD: this is done in Fig. 9 in the form of 
a scale invariant plot. We omit the errorbars for the sake of clarity, and just remind the 
reader that the errors - both for the two and the four fiavor model - are reasonably small 
up to hb — 500. We see that the transition in four flavor QCD hes consistently below that 
of the two flavor theory, and that this effect increases with increasing density. It looks as 
if the production of real fermions further favors the phase transition, which is indeed the 




(5.5) 



of ik/%r. 
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expected behaviour (see e.g. [36]). 



VI. CHIRAL CONDENSATE 

Taylor expansion and Fonricr decomposition are natural parametrisation for our ob- 
servables. In particular, the analysis of the phase diagram in the temperature-imaginary 
chemical potential plane suggests to use Fourier analysis for T < Tf.. Moreover, the fugacity 
expansion 

Z = ^c„e3'^"^* (6.1) 

hints that this would be easier in the cold phase, at stronger coupling, where a few coeffi- 
cients might suffice. At high temperature, in the weak coupling regime, on the other hand, 
perturbation theory might serve as a guidance, suggesting that the first few terms of the 
Taylor expansion might be adequate in a wider range of chemical potentials. 

As the chiral condensate is an even function of the chemical potential, its Fourier de- 
composition reads: 

(^^) = E «F cos{nNtN,i2i) (6.2) 

n 

which is easily continued to real chemical potential 

(ViV^) = E «F cosh{nNtN,fiR) (6.3) 

n 

In our Fourier analysis of the chiral condensate results we limit ourselves to n = 0, 1, 2 
and we assess the vahdity of the fits via both the value of the and the stability of ap and 
a]r given by one and two cosine fits. We summarize the results of the Fourier analysis of the 
chiral condensate results in Table IV. 

The Fourier analysis turns out to be satisfactory aX (3 = 5.01 and (3 = 5.03 where -as 
discussed in Section III and shown in Fig. 2 above- the chiral condensate is a continuous 
function of (3. One cosine fit is actually enough to describe our data at /? = 5.01 and 
P — 5.03, with the current statistical accuracy: adding a term cos{2NtNciJ,) — cos(24//) 
in the expansion does not modify the value of the coefficients and ap, and does not 
particularly improve the x^- The term cos(24yu) is anyway needed in order to asses the 
errors on the analytic continuation to real /i^, as we will discuss below. 

At /3 = 5.10 (in the fits we discarded of course the points corresponding to the RW 
discontinuity) we know (see again Section III and Fig. 2) that the periodicity is no longer 
smooth. Indeed, the values of the first two Fourier coefficients depend on the type of the fit 
(two or three parameters). 

Next, we have considered polynomial fits. In these fits we exploited the symmetries 
and the periodicity of the model to improve the statistical accuracy: all of the results were 
traslated, or symmetrised, to the first half period. For a quick comparison with the Fourier 
analysis we use a polynomial fit of the form 

(^^) = a^ - 72a^i/ + 864af /x^ (6.4) 
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When one cosine fit is adequate, the Taylor expansion would give — o}p-\-h\, and = af. 
At /3 = 5.01 and j3 = 5.03 this is indeed the case, within our largish errors. For = 5.10 a 
second order Taylor expansion is adequate at < 0.2, while a fourth order term does not 
substantially improve the behavior. As discussed above, the quality of the polynomial fits 
should improve at higher temperature, closer to the perturbative regime, where a second 
order polynomial should become exact. Results for the polynomial fits are collected in 
Table V. 

We can now continue the results of the Fourier analysis to real chemical potential. Fig. 
10 shows the behavior of the chiral condensate as a function of real chemical potential for 
/3 = 5.01 and (3 — 5.03, using 6.3 for both one and two cosine fits. We show the errorband 
induced by the fits with a cos(24/i) term in a large chemical potential range. The symbols 
(triangles and squares) are instead plotted only in the broken phase, i.e. for ji < ndP) (see 
next Section for more on this point). 

We have checked that the analytic continuation to real chemical potential of the results 
of the polynomial fits: 

= a^ + 72a^n^ + 864af / (6.5) 

is affected by comparable, if not smaller, systematic errors as the ones observed with the 
Fourier parametrization. 

As a final comment on the polynomial fit results, we consider the temperature dependence 
of a^. The Maxwell relation 

dm d/j, 

shows that a|' is the derivative with respect to the quark mass of the quark number suscepti- 
bility (by taking the derivative of cither sides of the relation above). These results sTiggcst 
that such derivatives -in lattice units- do not change much with temperature, in contrast 
with the quark number susceptibility itself: the results for the number density (we postpone 
a discussion of thermodynamics to another publication and we just quote the central values) 
are: 

Jo = 0.055/i+ 1.30;u3 (6.7) 

Jo = 0.073/i + 1.4 V (6.8) 
Jo = 0.455/x + 1.3 V (6.9) 

at P — 5.01,5.03,5.10 from top to bottom. The coefficient of the first term is the quark 
number susceptibility in lattice units, which increases rapidly in the plasma phase as it 
should. 



VII. NATURE OF THE CHIRAL PHASE TRANSTITION 

We have discussed in Section II the interrelation between chiral condensate and Polyakov 
loop at least up to = .2. Obviously, the observed correlation holds true at real chemical 
potential as well: we can consider the difference between the critical coupling for the chiral 
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condensate and the Polyakov loop: —l3^{n). Our results suggest that = 

at imaginary in a nonzero interval. It should then remain zero also at real chemical 
potential. We can then conclude that, within the present accuracy, the nearly coincidence of 
chiral and deconfining transition persists at nonzero baryon density in a four flavor model. 

As for the order of the phase transition, let us consider again Fig. 10 where, together 
with the errorband in a larger /j, interval, we have plotted , as triangles and squares, the 
results from one and two cosine fits, for fi < fidP), inferred from the results of Section III 
above. The behavior shown in Fig. 10 is then consistent with a first order transition. 

VIII. SUMMARY 

We have studied the phase diagram of four flavor QCD in the imaginary chemical po- 
tential - temperature plane. We have measured the location of the endpoint of the Roberge 
Weiss line and inferred that the RW line ends on the chiral critical line. 

We have continued our results to real chemical potential. The critical line in lattice units 
reads 

Pc = 5.037(3) - 0.96(10)/^^ + 0(/x^) 
which, in physical units, is well described by 

Te = yJ{T^ - 0.02 V) 

or, equivalently, by 

^ ^ / 0.021 .\ 

over a large interval of chemical potentials. It is of some interest to notice that this fea- 
ture is observed also in simple models, where we have found that the critical line is well 
approximated by a second order polinomial up to fiquark — T. 

We have found a good agreement with the results by Fodor and Katz [8] up to yu, ~ 400 
MeV. Our results are still compatible with theirs, within largish errors, up to ~ 500. 
Beyond that /i value our results do not have much statistical significance and a comparison 
is no longer meaningful. We emphasize that the main source of uncertainty on our results 
is statistical, and it is associated to a poor knowledge of the fourth order coefficient, which 
appears to be very small. 

We have studied the flavor dependence of the results by comparing our flndings for the 
critical line of four flavor QCD with those obtained by De Forcrand and Philipsen [16] in 
the two flavor model, and found it to be that expected on physical grounds. 

We have studied the character of the chiral transition. We have found indications that it 
remains correlated with the transition associated with the Polyakov loop. The dependence 
of the chiral condensate on the chemical potential is consistent with a flrst order transition. 
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TABLES 



TABLE L Locations of the chiral critical line 













Pc 










0.00 




5.0400(30) 










0.10 




5.0470(15) 










0.15 




5.0540(15) 










0.173(3) 




5.0650 










0.20 




5.0765(20) 










0.222(3) 




5.0850 










0.2617994 




5.0970(20) 








TABLE II. Results of least 


quadratic fits to the critical line 




Discarded 




a 


b 


xVd.o.f. 






0.00 




5.035(2) 


1.014(74) 


2.67 






0.10 




5.034(3) 


1.024(90) 


2.98 






0.15 




5.038(1) 


0.954(30) 


0.49 






0.175 




5.036(3) 


0.984(77) 


3.49 






0.20 




5.036(3) 


0.977(77) 


3.33 






0.22 




5.037(3) 


0.926(117) 


3.09 








TABLE III. 


Results of the least 


squares quadratic fits to the critical line 








range 


a 


b 


SSR 




pl 




[0:0.17] 


5.039(2) 


0.79(12) 


1.4C-05 


8 e-06 


p2 




[0:0.20] 


5.038(2) 


0.90(10) 


2.8C-05 


1.0 e-05 


p3 




[0:0.22] 


5.038(2) 


0.94(7) 


3.1e-05 


1.2 e-05 
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TABLE IV. Fourier coefficients for the chiral condensate. The constant and the first coefficeint 



are satisfactory determined for T < Te 



p 


n 


1 


9 

ap 


9/7 P 


5.03 


0.931(1) 


-0.0251(22) 


0(0) 


1.31 


5.03 


0.930(2) 


-0.0230(26) 


-0.0034(26) 


1.15 


5.01 


0.974(1) 


-0.0106(9) 


0(0) 


0.87 


5.01 


0.974 (1) 


-0.0107 (10) 


-0.0003 (12) 


0.98 


5.10 


0.396 (1) 


-0.024 (1) 


0(0) 


1.22294 


5.10 


0.404(4) 


-0.38(6) 


0.007(3) 


0.96 



TABLE V. Coefficients of a polynomial fit for the chiral condensate f(x) = 



- 72af x2 + 864af x^. 



p 


ao 


af 


af 


xVd.o.f. 


5.03 


0.905(2) 


-0.0272(40) 


-0.0261(70) 


1.17 


5.01 


0.963(1) 


-0.011(3) 


-0.009(4) 


1.03 


5.10 


0.372 (1) 


-0.019 (1) 


0(0) 


1.03 
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FIG. 1. Average value of the Polyakov loop phase as a function of the imaginary chemical 
potential for different values of /3. The vertical dashed line corresponds to 6 = m/T = 7r/3. 
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FIG. 2. Average value of the chiral condensate as a function of the imaginary chemical potential 
for different values of /?. The vertical dashed lines correspond to 9 = fii/T = {2k + 1)tt/3. The 
continuous line in the lower picture is the result of a quadratic fit at small values of a/x/ obtained 
at /3 = 5.10. 
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FIG. 3. Imaginary part of the barion density as a function of fij for different values of (3 
(left-hand side), and as a function of /? at 6* = /x//r = (right-hand side). 



20 



0.9 I 0.20 



° o 2 S 
0.8 5 D ^ ^ i ^ 



0.5 . 



P 



a 



0.15 



S 0.7 

fl o 

g 0.10 ® 

O pis 

2 0.6 1 

JS * 



0.05 



0.4 I 0.00 
5.07 5.08 5.09 5.1 5.11 5.12 5.13 



FIG. 4. Chiral condensate (white circles) and absolute value of the Polyakov loop (black circles) 
as a function of /3 for f? = vr/S. The sharp changes of the two quantities coincide with the location 
of endpoint of the RW critical line. 
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FIG. 5. Sketch of the phase diagram in the /Lt/-/3 plane. The filled cirles represents direct 
determinations of the chiral critical line location from our simulations. The rest of the chiral line 
has been obtained by interpolation and by exploiting the symmetries of the partition function. 
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FIG. 6. The critical line of four flavor QCD, ruq = .05, in the fj?, P plane. The dotted lines 
show the range of the results of the quadratic fits described in the text. The triangles are our results 
for the critical coupling from Table I. The results by Fodor and Katz (diamonds) are plotted for 
comparison. 
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FIG. 7. 0(/i^) effects on the critical line of four flavor QCD. The solid line shows the results 
of a second order fit, the thin dotted line the results of a fourth order fit with first and second 
order coefficients constrained to their central values. The upper (lower) thick dotted line shows the 
results of a fourth order fit with the second coefficent constrained to its best value plus (minus) its 
error. See text for details. 
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FIG. 8. The critical line of four flavor QCD, with ruq = .05. Result and error band are from 
P = 5.037(3) +0.96(10)/x^ converted to physical units. The result (central line) is well approximated 
by T = r,^(l-0.02lg. 
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FIG. 9. Scale invariant plot for the QCD critical line at high temperature, from imaginary 
chemical potential. Results for two flavors of staggered fermions (solid line) are from De Forcrand 
and Philipsen. The four flavor critical line (dotted) lies consistently below the two flavor one, the 
effect growing larger at larger density. 
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FIG. 10. The analytic continuation to real chemical potential of the chiral condensate: the large 
bands are the errorband associated with the second Fourier coefficient. The symbols - triangles 
and squares - are printed only in the broken phase, for one (smaller errors) and two (larger errors) 
Fourier coefficient results. The critical chemical potential - where the symbols stops- has been 
inferred from Fig. 8 and related discussions. The results are consistent with a first order chiral 
transition 



27 



